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ABSTRACT 

The neutrino-antineutrino annihilation into electron-positron pairs near the 
| surface of compact general relativistic stars could play an important role in 

OO . supernova explosions, neutron star collapse, or for close neutron star binaries 

near their last stable orbit. General relativistic effects increase the energy de- 
position rates due to the annihilation process. We investigate the deposition of 

Q\ ' energy and momentum due to the annihilations of neutrinos and antineutrinos 

O ' 

in the equatorial plane of the rapidly rotating neutron and quark stars, respec- 
tively. We analyze the influence of general relativistic effects, and we obtain 
' the general relativistic corrections to the energy and momentum deposition 



rates for arbitrary stationary and axisymmetric space-times. We obtain the 
energy and momentum deposition rates for several classes of rapidly rotating 
neutron stars, described by different equations of state of the neutron matter, 
and for quark stars, described by the MIT bag model equation of state and 
in the CFL (Color-Flavor-Locked) phase, respectively. Compared to the New- 
tonian calculations, rotation and general relativistic effects increase the total 
annihilation rate measured by an observer at infinity. The differences in the 
equations of state for neutron and quark matter also have important effects on 
the spatial distribution of the energy deposition rate by neutrino-antineutrino 
annihilation. 



Key words: neutrinos: dense matter - equation of state: stars: rotation: 
relativity. 
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1 INTRODUCTION 

Since the pioneering works of Cooperstein et al. (1986), Cooperstein et al. (1987), and 
Goodman et al. (1987), the energy deposition rate from the v + v — > e + + e~ neutrino 
annihilation reaction has been intensively studied. This reaction is of considerable impor- 
tance for Type II supernova dynamics, neutron star collapse, or for close neutron star bi- 
naries near their last stable orbit. Neutrino-antineutrino annihilation into electrons and 
positrons can deposit more than 10 51 ergs above the neutrino-sphere of a type II supernova 
(IGoodman et al. 19871 ). This energy deposition, together with neutrino-baryon capture, sig- 



nificantly in creases the neutrino 



l eafing in th e envelope via the so-called delayed shock 



mechanism ( Bethe fc Wilson 1985 



Bethe 19901 ) . For large r the energy deposition rate is 



proportional to r -8 , where r is the distance from the center of the neutrino-sphere. The ini- 
tial estimations of the neutrino annihilation reaction efficiencies were based on a Newtonian 
approach, by assuming that 2GM/c 2 R << 1, where M is the gravitational mass of the star, 
and R is the distance scale. However, for a full understanding of the effects of the neutrino 
annihi lation in strong gravitation al fields, general relativistic effects must be taken into ac- 
count (jSalmonson &: Wilson 19991 ) . For a static neutron star, by adopting for the description 
of the gravitational field the Schwarzschild metric, the efficiency of the v + v — > e + + e~ 
process is enhanced over the Newtonian values up to a factor of more than 4 in the regime 
applicable to Type II supem ovae, and by up to a factor of 30 for collapsing neutron stars 
(jSalmonson fc Wilson 19991 ) . The neutrino pair annihilation rate into electron pairs between 
two neutron stars in a binary system was calculated by Salmonson & Wilson (2001). A closed 
formula for the energy deposition rate at any point between the stars was obtained, where 
each neutrino of a pair derives from each star, and this result was compared with that in 
which all neutrinos derive from a single neutron star. An approximate generalization of this 
formula was also given to include the relativistic effects of gravity. The interstar neutrino 
annihilation is a significant contributor to the energy deposition between heated neutron 
star binaries. 

The neutrino-antineutrino annihilation into electrons and positrons is a n important can- 
didate to explain the energy source of the gamma ray bursts (GRBs) (jPaczvns ki 1990; 



Meszaros fc Rees 1992 



Ruffert fc Janka 1998 



Ruffert fc Janka 1999 



Asano fc Iwamoto 20021 ) 



* E-mail: zkovacs@mpifr-bonn.mpg.dc 
f E-mail: hrspksc@hkucc.hku.hk 
| E-mail: harko@hkucc.hku.hk 



Electron-positron energy deposition near neutron and quark stars 3 

The semi-analytical study of the gravitational effects on neutrino pair annihilation near the 
neutrinosphere and around the thin accretion disk were considered in Asano & Fukuyama 
(2000), by assuming that the accretion disk is isothermal, and that the gravitational field 
is dominated by the Schwarzschild black hole. General relativistic effects were studied only 
near the rotation axis. The energy deposition rate is enhanced by the effect of orbital bend- 
ing toward the center. However, the effects of the redshift and gravitational trapping of the 
deposited energy reduce the effective energy of the gamma-ray burst's source. Although each 
effect is substantial, the effects partly cancel one another. As a result, the gravitational effects 
do not substantially c hange the energy depositio n rate for either the spherically symmetric 
case or the disk case (jAsano &: Fukuyama 20001 ). Using idealized models of the accretion 
disk, Asano & Fukuyama (2001) investigated the relativistic effects on the energy deposition 
rate via neutrino pair annihilation near the rotation axis of a Kerr black hole, by assuming 
that the neutrinos are emitted from the accretion disk. The bending of neutrino trajectories 
and the redshift due to the disk rotation and gravitation were also taken into consideration. 
The Kerr parameter, a, affects not only behavior of the neutrinos, but also the inner radius 
of the accretion disk. When the deposition energy is mainly contributed by the neutrinos 
coming from the central part, the redshift effect becomes dominant as a becomes large, 
and the energy deposition rate is reduced compared with that neglecting the relativistic 
effects. On the other hand, for a small a, the bending effect becomes dominant and makes 
the energy increase by fact or of 2, compared with that which neglects the relativistic effects 
(jAsano fc Fukuyama 200ll ). 

The effect of the inclusion of the slow rotation of the star in the general relativistic treat- 
ment of the neutrino-antineutrino annihilation into electron positron pairs was considered 
in Prasanna & Goswami (2002). It was shown that the inclusion of the rotation results in a 
reduction in the heating rate, as compared to the no rotation case. The energy-momentum 
deposition rate (MDR) from the v — v collisions above a rotating black hole/thin accre- 
tion disk system was calculated by Miller et al. (2003), by imaging the accretion disk at a 
specified observer using the full geodesic equations, and calculating the cumulative MDR 
from the scattering of all pairs of neutrinos and antineutrinos arriving at the observer. The 
dominant contribution to the MDR comes from near the surface of the disk with a tilt of 
approximately 7r/4 in the direction of the disk's rotation. The MDR at large radii is directed 
outward in a conic section centered around the symmetry axis and is larger by a factor of 10- 
20 than the on-axis values. There is also a linear dependence of the MDR on the black hole 
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angular momentum. The deposition of energy and momentum, due to the annihilation of 
neutrinos v and antineutrinos v in the vicinity of steady, axisymmetric accretion tori around 
stellar-mass black holes was investigated in Birkl et al. (2007). The influence of general rel- 
ativistic effects were analyzed in combination with different neutrinosphere properties and 
spatial distribution of the energy deposition rate. Assuming axial symmetry, the annihilation 
rate 4- vector was numerically computed. The local neutrino distribution was constructed by 
ray-tracing neutrino trajectories in a Kerr space-time, using null geodesies. Different shapes 
of the neutrinospheres, spheres, thin disks, and thick accretion tori were studied, whose 
structure ranges from idealized tori to equilibrium non-selfgravitating matter distributions. 
Compared to Newtonian calculations, general relativistic effects increase the total annihila- 
tion rate measured by an observer at infinity by a factor of two when the neutrinosphere is 
a thin disk, but the increase is only 25% for toroidal and spherical neutrinospheres. Thin 
disk models yield the highest energy deposition rates for neutrino-antineutrino annihilation, 
and spherical neutrinospheres the lowest ones, independently of whether general relativistic 
effects are included. General relativity and rotation cause important differences in the spa- 
tial distribution o f the energy deposition rate by neutrino v and antineutrino ^-annihilation 
( IBirkl et al. 20071 ). The study of the structure of neutron star disks based on the two-region 
(i.e., inner and outer) disk scenario was performed by Zhang & Dai (2009), who calculated 
the neutrino annihilation luminosity from the disk in various cases. The effects of the vis- 
cosity parameter a, energy parameter e (measuring the neutrino cooling efficiency of the 
inner disk), and outflow strength on the structure of the entire disk, as well as the effect of 
emission from the neutron star surface boundary emission on the total neutrino annihilation 
rate were investigated. An outflow from the disk decreases the density and pressure, but in- 
creases the thickness of the disk. Moreover, compared with the black hole disk, the neutrino 
annihilation luminosity above the neutron star disk is higher, and the neutrino emission from 
the boundary layer could increase the neutrino annihilation luminosity by about one order 
of magnitude higher than the disk without boundary emission. The neutron star disk with 
the advection-dominated inner disk co uld produce the hig hest neutrino luminosity, while 
the disk with an outflow has the lowest (jZhang fc Dai 20091 ). A detailed general relativistic 
calculation of the neutrino path for a general metric describing a rotating star was studied 
in Mallick & Majumder (2009). The neutrino path was calculated along the equatorial and 
polar plane. The expression for the minimum photosphere radius was obtained and matched 
with the Schwarzschild limit. The minimum photosphere radius was calculated for stars with 
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two different equations of state, each rotating with two different velocities. The results show 
that the minimum photosphere radius for the hadronic star is much greater than for the 
quark star, and that the minimum photosphere radius increases as the rotational velocity 
of the star decreases. The minimum photosphere radius along the polar plane is larger than 
that along the equatorial plane. The estimate of the energy deposition rate for neutrino pair 
annihilation for the neutrinos coming from the equatorial plane of a rotating neutron star 
was calculated along the rotation axis in Bhattacharyya et al. (2009), by using the Cook- 
Shapiro- Teukolsky metric. The neutrino trajectories, and hence the neutrino emitted from 
the disk, are affected by the redshift due to the disk rotation and gravitation. The energy 
deposition rate is very sensitive to the value of the temperature, and its variation along the 
disk. The rotation of the star has a negative effect on the energy deposition rate, it decreases 
with increase in rotational velocity. The standard model for Type II supernovae explosions, 
confirmed by the detection of neutrinos emitted during the supernova explosion, predict s 
the formation of a compact object, usually assumed to be a neutron star (IChan et al. 20091 ) . 
However, the newly formed neutron star at the center of SN 1987A may undergo a phase 
transition after the neutrino trapping timescale (~ 10 s). Consequently the compact rem- 
nant of SN 1987A may be a strange quark sta r, which has a softer equation of state than 
that of neutron star matter (IChan et al. 20091 ). Such a phase transition can induce stellar 
collapse and result in large amplitude stellar oscillations. A three-dimensional Newtonian 
hydrodynamic code was used to study the time evolution of the temperature and density at 
the neutrinosphere. Extremely intense pulsating neutrino fluxes, with submillisecond period 
and with neutrino energy (greater than 30 MeV), can be emitted because the oscillations 
of the temperature and density are out of phase almost 180°. The dynamical evolution of 
a phase-transition-induced collapse of a neutron star to a hybrid star, which consists of a 
mixture of hadronic matter and strange quark matter, was studied in Cheng et al. (2009). It 
was found that both the temperature and the density at the neutrinosphere are oscillating 
with acoustic frequency. Consequently, extremely intense, pulsating neutrino/antineutrino 
fluxes will be emitted periodically. Since the energy and density of neutrinos at the peaks 
of the pulsating fluxes are much higher than the non-oscillating case, the electron/positron 
pair creation rate can be enhanced dramatically. Some mass layers on the stellar surface can 
be ejected, by absorbing energy of neutrinos and pairs. These mass ejecta can be further ac- 
celerated to relativistic speeds by absorbing electron/positron pairs, created by the neutrino 
and antineutrino annihilation outside the stellar surface. 
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It is the purpose of the present paper to consider a comparative systematic study of 
the neutrino-antineutrino annihilation process around rapidly rotating neutron and strange 
stars, respectively, and to obtain the basic physical parameters characterizing this process 
(the electron-positron energy deposition rate per unit volume and unit time, and the total 
emitted power, respectively), by taking into account the full general relativistic corrections. 
In order to obtain the electron-positron energy deposition rate for various types of neutron 
and quark stars we generalize the relativistic description of the neutrino-antineutrino anni- 
hilation process to the case of arbitrary stationary and axisymmetric geometries. To com- 
pute the electron-positron energy deposition rate, the metric outside the rotating general 
relativistic stars must be determined. In the present study we study the equilibrium con- 
fig urations of the rotating neutro n and quark stars by using the RNS code, as introduced 



in (IStergioulas fc Friedman 19951 ). and discussed in detail in (IStergiou 



was used for the study of different models of rotating ne utron stars in (Nozawa et al. 19981 ) 



as 2003). This code 



and for the study of the rapidly rotating strange stars ( (IStergioulas et al. 19991 ) ). The soft- 
ware provides the metric potentials for various types of compact rotating general relativistic 
objects, which can be used to obtain the electron-positron energy deposition rate in the 
equatorial plane of rapidly rotating neutron and quark stars. 

The present paper is organized as follows. In Section 2 we present the basic formalism for 
the calculation of the electron-positron energy deposition rate from neutrino-antineutrino 
annihilation. The general relativistic corrections to this process are obtained in Section 3. 
The equations of state of dense neutron and quark matter used in the present study are 
presented in Section 4. In Section 5 we obtain the electron-positron energy deposition rates 
in the equatorial plane of the considered classes of neutron and quark stars. We discuss and 
conclude our results in Section 6. 



2 NEUTRINO PAIR ANNIHILATION 



A considerable amount of energy can be released by the neutrino pair annihilation process 
in the regions close to the so called neutrino-sphe re, with radius R v , at which the mean free 
path of the neutrino is equal to the radius itself (IBethe fc Wilson 19851 ). 
The energy deposition rate per unit volume is given by 



q(r) 



fuiPu, r)fi>(p D ,r){a\v v - v p \e u e p } 



£j ^d 3 Pu d 3 p, 



(1) 



for any point r > R u , where f v and fp are the number densities in the momentum spaces 
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with the momenta p v and p_, and v v , e v , v D and e D are the 3- velocities an d the energy of the 
colliding neutrino-antineutrino pairs, respectively (IGoodman et al. 19871 ). The cross-section 
of the collision is denoted by a. 

By applying the decompositions p u = e u Q. v and d z p u = elde u dfl u , with the solid angle 
vector Q, v pointing in the direction of p u , and with the assumption that the neutrino-sphere 
emits particles isotropically, the integral in Eq. ([I]) can be separated into an energy integral 
and an angular part. After evaluating the energy integral for fermions, for a spherically 
symmetric geometry the energy deposition rate per unit volume can be represented as 

«K0 = 7 ^ffi (5 W e(r) cx T 9 (r)0(r), (2) 

where T is the neutrino temperature, £ is the Riemann function, D = 1±4 sin 2 #iy+8 sin 4 9\y, 
with sin 2 8\y = 0.23, and G 2 F = 5.29 x 1CT 44 cm 2 MeV -2 , respectively, and with the angular 
part of the energy deposition rate 0(r) given by 

0(r) = JJ{1-Sl v - fl p ) 2 dfl u dfl p . (3) 

The radial momentum density p transported into the e + e~ plasma from the colliding vv 
pairs can be written as (jsalmonson fc Wilson 1999 ) 

p = U p (r)Q(r), (4) 
where 

%{r) = -j J \l - n v ■ fi p ) 2 [r • (a, + ft p )]dfl u dfl p , (5) 

with the unit vector r normal to the stellar surface. 

Since the integrals in Eqs. (jSJ) and (jSJ) depend only on the radial coordinate r, by virtue 

of the symmetry, the angular part G(r) and the function $ p (r) can be given by the analytic 

formulae 

2vr 2 

0(r) = — (l-x) 4 (x 2 + 4x + 5) (6) 

o 

and 



7T 



2 



$ p ( r ) = _(1 - x ) 4 (8 + 17x + 12x 2 + 3x 6 ) , (7) 
where x = cos9 mnx , and is the maximal angle between r and p v ( IGoodman et al. 1987 ; 



Salmonson fc Wilson 19991 ) . In order to determine the quantity x, one need to use the equa- 



tions of motion of the neutrinos radiated by the stellar matter, and propagating outside 
the neutrino-sphere. The equations of motion in the Newtonian case, or the geodesic equa- 
tions in the general relativistic case fully determine x as a function of the radial coordinate. 
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Therefore, the properties of the gravitational potential produce an imprint on the energy 
deposition rate calculated in the region close to the neutrino-sphere. Salmonson & Wilson 
(1999) extended the calculations of Goodman et al. (1987) to the case of the Schwarzschild 
geometry, and compared their results to those of the Newtonian case. 

Although rotating configurations of neutron stars break the spherical symmetry of the 
space-time, one can still carry out an analysis similar to the static case if the study is 
restricted to the investigation of the annihilation of neutrino and antineutrino pairs propa- 
gating in the equatorial plane of the rotating star. By eliminating the angular dependence 
from the equations, a formalism similar to the spherically symmetric case can be used to 
calculate the energy deposition rate in the equatorial plane. The obtained result is not equal 
to the total deposition rate of the high energy electron positron pairs created in the annihi- 
lation process. However, this quantity can still be applied in the comparison of the neutrino 
and antineutrino annihilation energy deposition rate for different models of rotating neutron 
and quark stars, or in studying the general effects of the rotation of the stellar object on 
this process. 



3 GENERAL RELATIVISTIC EFFECTS ON THE 

ELECTRON-POSITRON ENERGY DEPOSITION RATE 

The metric of a stationary and axisymmetric geometry is given in the general form by 

ds 2 = g tt dt 2 + 2g t( pdtd(j) + g rr dr 2 + geedO 2 + g^dcj) 2 . (8) 

For this metric the null-geodesics equations in the equatorial plane 6 = n/2 are 

^ = 9h E + 9i4 L ( 9 ) 

9% ~ 9tt9<P4> ' 
• g<p t E + g tt L 

= 2 > ( 10 ) 

9t<t> ~ 9ttg<t>4> 

. 2 g H E 2 + 2g t4> EL + g tt L 2 
g rr r — 2 , ll ) 

9t<f> ~ 9ttg<t>4> 

where E is the energy, and L is the angular momentum of the particles propagating along 
the null-geodesics. Then the parametric equation for dr/d(f) can be written as 

/ dr \ 2 _t„2 „ n 9 w + 2 9t<pb + g tt b 2 

{9t<t> + 9ttby 

where the impact parameter is defined as b = L/E. 



9rr(7l) ={9%-9u9h) 7——TT 2 , (12) 
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For the metric given by Eq. ([81) , the locally non- rotating frame (LNRF) (in which the 
wordlines of the freely falling observers are r = constant, 9 = constant and <j)—ujt = constant, 
respectively, with u = —gt^/g^) has the basis of one- forms (jBardeen et al. 1972T ) 

Jt) _ 



-9a 



-1,0,0,0), 



(r) 
(«) 
(0) 



\ 9t<t> ~ 9tt9<p<t> 

(0, y/fo,0,0), 

(0,0,^,0), 

Since the velocity measured in the LNRF is given by = eJV, the angle 6* r between 
the particle trajectory and the tangent vector to the circular orbit with the radial coordinate 
r can be written as 

e^v r Jlfrr dr 



tan 6L 



V W e W^ + e W y /m[ 1 + 9t4>/(9H v *)]d<l>' 
From this expression we obtain for dr/d<j) the equation 

9t<t> 



dr 



1 + 



tan r 



a<P V V 9u v ) 
By inserting Eqs. ([9]) and (TT01) into the definition of f*, we obtain for the expression 

* g t<t> E + g tt L 



(13) 



u 1 



which can be substituted into Eq. ([TBI) to give 



1 - 



tan 6L 



a<p y # rr \ g^ g tl pE + g tt L j 
Then the derivative dr/dip can be eliminated from the parametric equation (1121) and we 
obtain 



9t<t> 9<j><t> + 9t4>b 
9t</> + 9ttb , 



tan 2 9 r = (gL - g tt 9<w) 



+ 2g t(f) b + g tt b 2 
(,9<j>t + 9ttb) 2 



This result gives a second order algebraic equation for b, 



[(9u9h> ~ 9^) sec2 °r + 9^b 2 + 2g 1 



+ 



0. 



(14) 



which can be solved to give the impact parameter b as 
b± 



= • (15) 

9t<t> ± \j9t4, ~ gttg^sec e r 

A particular system of coordinates that is used in the study of the general-relativistic 

rotating configurations is the q uasi- isotropic coordinate system (t,r,6,4>), in w hich the line 

element can be represented as ( jStergioulas fc Friedman 19951 ; IStergioulas 20031 ) 
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ds 2 = -e^ +p dt 2 + e 2a (df 2 + f 2 d9 2 ) + e^f 2 sin 2 9 {dxj> - Qdtf 



(16) 



where 7, p, a and the angular velocity of the stellar fluid relative to the local inertial frame 
us are all functions of the quasi-isotropic radial coordinate f and of the polar angle 9. 

If for neutron stars the metric (jSJ) is given in an isotropic coordinate system in the form 
(j!6p . then the second order algebraic equation (j!4p for h can be written in terms of the 
metric functions p and Q as 

(sec 2 9 r + e- 2p u) 2 f 2 sin 2 9) b 2 - 2e~ 2p Qf 2 sin 2 9b + e~ 2p f 2 sin 2 = 0. 



For the impact parameter, corresponding to 9 r = 0, we obtain (ICadeu et al. 20071 ) 
— e~ p fsm9 e~ p fsm9 



—e~ p uor sin 9 ± sec # r sec 9 r ± e~ p tur sin 
From the parametric equation Eq. ( 1121) we obtain the deflection angle of the particle 
trajectory for a given 6 as 

y/gVr(gt<f> + 9ttb)dr 



fobs 



(17) 



- 9u9h>)(9h> + 2 9t<t>b + gttb 2 ) 
In this equation measures the change in the angle between the source and the ob- 
server, for a photon emitted at the radial coordinate r em , and observed at the radial coor- 
dinate r b s . This equation can also be given in the equatorial plane in terms of the metri c 
functions p and Co, respectively, appearing in the line element Eq. fflBT) (ICadeu et al. 20051 ) 



A0 



f ° bs c a-(7+p)/2 "{1 



u)b) + be 2p /f 



-.dr. 



- cub) 2 - b 2 e 2p /f 2 

In the equatorial plane of black holes the photon radius is defined as the innermost 
bounda ry of circular orbits below which massless particles with 9 r = are gravitationally 
bound (IBardeen et al. 19721 ). For static black holes the photon radius is 3M, and for rotating 
black holes it reduces to M, as the spin parameter a* = J/M 2 of the black hole approaches 
unity. This orbit may exist for ultra-compact stars as well. For static stars with a stellar 
radius less than 3M there is always such a "photon sphere", whereas, depending on the 
geometry of the s pace-time, the rotating star s can also have a photon radius at both lower 
and higher radii (jMallick fc Majumder 20091 ) . On the other hand, very massive rotating 
quark stars in the Color-Flavor-Locked phase can reach masses higher than the equilibrium 
limit for static stars (~ 3M ), of t he same order as th e stellar mass black holes, and thus 
they can also have a photon radius (IKovacs et al. 20091 ). 
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If the integral in Eq. (|T7|) diverges to infinity (or dr/dip tends to zero), the null particles 

are rotating around the central objects in circular orbits. In this case 9 r = 0, and the 

algebraic equation Eq. (HM for b reduces to + 2g t( j ) b + g u b 2 = 0, with the solution 

6± = -9u = ± e^r (lg) 

g t4> ± yfifo - g tt g H 1 ± e "Of 

For any value of b satisfying Eq. (|T8j) . the integral (fT7|) is divergent, and the null particles 
have circular orbits. In Eq. (fT8l) the impact parameter in the equatorial plane is given as a 
function of the radial coordinate only. 

In the case of the ultracompact static stars with radii R e less than 3M the (local) 
maximum of the function b(r) is located at r = 3M, providing the photon radius. If R e > 
3M, then b is a monotonically decreasing function, without a local maximum. For rotating 
compact stars, the function b(r) provides the same criterion for the existence of the potential 
barrier: for massless particles the equatorial orbit where b(r) attains its local maximum 
defines the innermost boundary of circular orbits. Even if this value is less than the equatorial 
radius R e , neutrinos are still free to propagate along orbits lying on the photon radius. 

By assuming that the mean free path of the neutrinos is equal to or less than the photon 
radius, Salmonson & Wilson (1999) identified the photon and the neutrino spheres with each 
other. Accordingly, we will also consider the orbit at the photon radius as the minimal radius 
where the annihilation process should still be taken into account, provided it is outside the 
star. The contribution of the electron-positron pairs formed inside the star to the deposition 
rate is neglected, because of their complicated interactions with the neutron and quark 
matter. 

Since the impact parameter measured at infinity is constant along the trajectory of any 
null particle, the neutrinos propagating from a point at the photon sphere radius R (or the 
stellar surface R e ) with the angle Or = tan -1 (R) / 'v^ (R) , will reach another point with 
radial coordinate r with the angle 



COSV r = F , q . (19) 

9M>{r) [gt^R) + \Jg?4>(R) - gu{R)g^{R) sec0 fl J - g^{R)g t <t>{r) 

Eq. f)19p allows to express x in the analytic expression of the angular part 0(r) of the energy, 
given by Eq. (jBJ), as 

2/ x , 9U( R ) fc( r ) ~ 9tt{r)g^{r)] 
x {r) = l- -2. (20) 



{g^( r ) [gt<t>{R) + ^/gt4 R ) - 9tt(R)g<M>(R)\ - g<M>(R)9t*{r)} 

For the line element given by Eq. ([TBI) . Eq. (120]) has the form 
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x 2 (r) = 1 



[fl + [u(f)-u(R)]Re^j 
Thus the function B(r) can be represented in the general form 



(21) 



6(r) 



2tt 2 



qU r ) 


9t^( r ) - gtt{r)g H {r)_ 




{g^( r ) 


g t <t>{R) + ^/gU R ) - 9u(R)g^(R)_ 


1 2 



gl4 R ) fe( r ) _ 9tt{r)g^{r) 



\g<t><t>(r) [gt<t>(R) + \jg^{R) - gtt(R)g^(R)] - g^(R)gt<p(r)} 



+4 



1 - 



g^iR) gttir) - gtt(r)g^(r) 



\ {g^(r) [gt^(R) + ^Jgt4 R ) - 9tt{R)g<M>( R )\ ~ gu{ R )9t4>( r )\ 

In quasi-isotropic coordinates we have 



+ 6 



• (22) 



6(f) 



2tt 2 



x 



1 - 


-{ 


\ 





R 



o 7 (f)- 7 (R) 



R 



r i _( 

e 7(f)- 7 (H) 



r i + 



u(f) - lo(R) 



uj(r)-uj(R)) Re^ R ) 

2 

+ 4. 



Re^ R ) 



\ 



1 



R 



, 7 (f)-7(R) 



r 1 + 



) 



+ 6 . 
(23) 



The neutrino temperature at the radius r can be expressed in terms of the temperature of 
the neutrino stream at the neutrino-sphere radius R, by taking into account the gravitational 
redshift. The redshift formula for T is the same as for the photon energy, 



T(r) 





'g^(R) - g u {R)g H {R) 




g^( r ) - gtt{r)g H (r) 



1/2 



T{R). 



(24) 



For the observed luminosity of the neutrino annihilation the redshift relation is given by 



g H {r -> oo) 


g%{R) - gtt(R)g^{R) _ 




g<t»t>( R ) 


g 2 t^ r oo) - gtt(r -> oo)g H {r - 


-> oo) 



r ^2 



L(R) 



-gu(R) 



L{R) , (25) 



_g<i>4>{ R ) 

since for isolated gravitating systems, such as rotating stars, the spacetime is asymptotically 
flat. Here the neutrino luminosity at the neutrino-sphere is 



L(R) = L U + L P = {A7iR 2 )^-acT 4 (R) 

lb 



(26) 



where a is the radiation constant. In this formula the curvature radius R is used to obtain 
the total area of the spherical surface through which the neutrino radiation is emitted. If 
we insert Eq. (1201) . describing the path-bending of the neutrinos, and the redshift formulae 
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Eqs. fl24|) - fl26j) into the decomposed expression Eq. ([2]) of q, we can calculate the effects of 
the gravitational potential on the deposition rate in the equatorial plane: 

{/ ^ 9/2 

^fc^ R-*\ (27) 
\J9u( R ) [9t^( r ) ~ 9tt(r)g^(r)\ J 

The proportionality factor, omitted from Eq. (I27"j) . is the same as the one in the New- 
tonian case, which will be used in the following to normalize the deposition rate for the 
general relativistic case. Eq. 0271) describes the energy deposition rate in e + e~ pairs from 
the neutrino-antineutrino annihilation process at radius r in the equatorial plane above the 
neutron or quark star photon sphere radius R, and with the neutrino luminosity observed 
at infinity L^. This relation can be also given in the quasi-isotropic coordinate system in 
terms of the metric functions of the line element Eq. ( jl~6l) . 

q(f) oc L 9 / 4 6(f)e 9 ^^ + ^^]/ 4 - 9 [^ f ) + ^ f )]/ 2 i?- 9 / 4 (/?). (28) 



4 EQUATIONS OF STATE AND STELLAR MODELS 

In order to obtain a consistent and realistic physical description of the rotating general 
relativistic neutron and quark stars, as a first step we have to adopt the equations of state 
for the dense neutron and quark matter, respectively. In the present study we consider the 

following equations of state for neutron and quark matter: 

1) Akmal-Pandharipande-Ravenhall (APR) EOS flAkmal et al. 19981 ) . EOS APR has 
been obtained by using the variational chain summation methods and the Argonne v\s two- 
nucleon interaction. Boost corrections to the two-nucleon interaction, which give the leading 
relativistic effect of order (v/c) 2 , as well as three-nucleon interactions, are also included in 
the nuclear Hamiltonian. The density range is from 2 x 10 14 g/cm 3 to 2.6 x 10 15 g/cm 3 . The 
maximum mass limit in the static case for this EOS is 2.20M^. We jo in this equation of 
state to the composit e BBP (e/c 2 > 4.3 x 10 n g/cm 3 ) faavm et al. 1971al) - BPS (10 4 g/cm 3 



< 4.3 x 10 n g/cm 3 ) flBavm et al. 1971bf ) - FMT (e/c 2 < 10 4 g/cm 3 ) flFevnman et al. 19491 ) 



equations of state, respectively. 



2) Douchin-Haensel (DH) EOS (jDouchin &: Haensel 20011 ) . EOS DH is an equation of 
state of the neutron star matter, describing both the neutron star crust and the liquid core. 
It is based on the effective nuclear interaction SLy of the Skyrme type, which is particularly 
suitable for the application to the calculation of the properties of very neutron rich matter. 
The structure of the crust, and its EOS, is calculated in the zero temperature approximation, 
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and under the assumption of the ground state composition. The EOS of the liquid core is 
calculated assuming (minimal) npefi composition. The density range is from 3.49 x 10 11 
g/cm 3 to 4.04 x 10 15 g/cm 3 . The minimum and maximum masses of the static neutron stars 
for this EOS are O.O94M and 2.O5M , respectively. 

3) Shen- Toki-Oyamatsu-Sumiyoshi (STOS) EOS (jShen et al. 19981 ). The STOS equation 
of state of nuclear matter is obtained by using the relativistic mean field theory with nonlin- 
ear a and uj terms in a wide density and temperature range, with various proton fractions. 
The EOS was specifically designed for the use of supernova simulation and for the neutron 
star calculations. The Thomas-Fermi approximation is used to describe inhomogeneous mat- 
ter, where heavy nuclei are formed together with free nucleon gas. We consider the STOS 
EOS for several temperatures, namely T = 0, T = 0.5 and T = 1.0 MeV, respectively. The 
temperature is mentioned for each STOS equation of state, so that, for example, STOS 
represents the STOS EOS for T = 0. For the proton fraction we chose the value Y p = 10~ 2 
in order to avoid the negative pressure regime for low baryon mass densities. 

4) Relativistic Mean Field (RMF) equations of state wi th isovector scalar mean fi eld cor- 
responding to the 5- meson- RMF soft and RMF stiff EOS fjKubis fc Kutschera 19971 ). While 
the 5-meson mean field vanishes in symmetric nuclear matter, it can influence properties of 
asymmetric nuclear matter in neutron stars. The Relativistic mean field contribution due to 
the 5-field to the nuclear symmetry energy is negative. The energy per particle of neutron 
matter is then larger at high densities than the one with no 5-field included. Also, the proton 
fraction of /3-stable matter increases. Splitting of proton and neutron effective masses due 
to the £-field can affect transport properties of neutron star matter. The equations of state 
can be parameterized by the coupling parameters C 2 = g 2 /m 2 , C 2 = g^,/m^, b = bjg\ and 
c = c/g*, where m a and m w are the masses of the respective mesons, and b and c are the 
coefficients in the potential energy U (er) of the cx-field. The soft RMF EOS is parameterized 
by Cl = 1.582 fm 2 , C 2 = 1.019 fm 2 , b = -0.7188 and c = 6.563, while the stiff RMF EOS 
is parameterized by C 2 = 11.25 fm 2 , C 2 = 6.483 fm 2 , b = 0.003825 and c = 3.5 x KT 6 , 
respectively. 



5) Baldo-Bombaci-Burgio (BBB) EOS flBaldo et al. 19971 ). The BBB EOS is an EOS 



for asymmetric nuclear matter, derived from the Brueckner-Bethe-Goldstone many-body 
theory with explicit three-body forces. Two EOS's are obtained, one corresponding to the 
Argonne AV14 (BBBAV14), and the other to the Paris two-body nuclear force (BBBParis), 
implemented by the Urbana model for the three-body force. The maximum static mass con- 



Electron-positron energy deposition near neutron and quark stars 15 

figurations are M max = 1.8M Q and M max = 1.94M when the AV14 and Paris interactions 
are used, respectively. The onset of direct Urea processes occurs at densities n ^ 0.65 fm~ 3 
for the AV14 potential and n ^ 0.54 fm~ 3 for the Paris potential. The comparison with 
other microscopic models for the EOS shows noticeable differences. The density range is 
from 1.35 x 10 14 g/cm 3 to 3.507 x 10 15 g/cm 3 . 



6) Bag model equation of s tate for quark matter (Q) EOS (jltoh 197CH ; 



Bodmer 1971 



Witten 1984 : 



Cheng et al. 19981 ). For the description of the quark matter we adopt first a 



simple phenomenological description, based on the MIT bag model equation of state, in 
which the pressure p is related to the energy density p by 



V 



(p - AB) c 2 



(29) 



where B is the difference between the energy density of the perturbative and non-perturbative 
QCD vacuum (the bag constant), with the value 4B = 4.2 x 10 14 g/cm 3 . 

7) It is generally agreed today that the color-flavor-locked (CFL) state is likely to be the 
ground s tate of matter, at least for asymptotic densities, and even if the quark masses ar e 



unequal (lAlford et al. 1999 



Rapp et al. 2000 



Horvath fc Lugones 2004 ; lAlford et al. 20071 ). 



Moreover, the equal number of flavors is enforced by symmetry, and electrons are absent, 
since the mixture is automatically neutral. By assuming that the mass m s of the s quark 
is not large as compared to the chemical potential /i, the thermodynamical pot ential of the 
quark matter in CFL phase can be approximated as (ILugones fc Horvath 2002l ) 



n 



CFL 



3/i 4 3m 



+ 



471-2 4^2 



12 In (m 8 /2/i) 
32vr 2 



A V + B, 



Ti 



(30) 



where A is the gap energy. With the use of this expression the pressure P of the quark 

matter in the CFL phase can be ob tained as an explicit function of the energy density e in 

the form (ILugones fc Horvath 2002l ) 
1 . . _ 2A 2 5 2 m 2 S 2 



P = -(e-4B) + 



where 
5 2 = - 



7T 



2tt 2 



a + \ a 2 + -7r 2 (e 



4 
9' 



B), 



(31) 



(32) 



and a = —m 2 /6 + 2A 2 /3. In the following the value of the gap energy A considered in each 
case will be also mentioned for the CFL equation of state, so that, for example, CFL200 
represents the CFL EOS with A = 200. For the bag constant B we adopt the value 4B = 
4.2 x 10 14 g/cm 3 , while for the mass of the strange quark we take the value m s = 150 MeV. 
The pressure-density relation is presented for the considered equations of state in Fig. [D 
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Figure 1. Pressure as a function of density (in a logarithmic scale) for the equations of state DH, RMF soft, RMF stiff, STOS 
0. STOS 0.5, STOS 1, BBBAV14, BBBParis, APR, Q, CFL150, and CFL300, respectively. 



To calculate the equilibrium configurations of the rotating neutron and quark stars with 
the EOS's presented here we use the RNS code, as introduced in Stergioulas & Friedman 
(1995), and discussed in details in Stergioulas (2003). This code was used for the study 
of different models of rotating neu tron stars (INozawa et al. 19981) . and for the study of 
the rapidly rotating strange stars (jstergioulas et al. 19991 ). The RNS code produces the 
metric functions in a quasi-spheroidal coordinate system, as functions of the parameter 
s — fj (f + f e ), where f e is the equatorial radius of the star, which we have converted into 
Schwarzschild-type coordinates r according to the equation r = f exp [(7 — p) /2]. 



5 ELECTRON-POSITRON ENERGY DEPOSITION RATE IN THE 
EQUATORIAL PLANE OF RAPIDLY ROTATING NEUTRON AND 
QUARK STARS 

To demonstrate the existence and the location of the photon radius for neutron and quark 
stars in the static and the rotating cases, respectively, in Fig. [2] we present b(r)/M as 
a function of the Schwarzschild coordinate radius (normalized to the mass of the star). 
For the static neutron star (modeled, for example, by the STOS 0.5 EOS), the exterior 
geometry is always the Schwarzschild one, and the curve b(r) connects to the curve of the 
static black hole at a radius somewhat lower than 3M. Then, the curve has a local minimum 
at 3M, similarly to the c ase of the function b(r), corresponding to the Schwarzschild case 
( iLandau fe Lifshitz 19981 ). Therefore, the star has a photon radius at 3M. This is not the 
case for the static quark star (with an EOS of the CFL 300 type), which is not so compact, 
and its curve b(r) reaches the Schwarzschild black hole curve at a somewhat higher radius 
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Figure 2. The impact parameter b/M as a function of the Schwarzschild coordinate r/M for static black holes, and static and 
rotating neutron and quark stars, respectively. 

than 3M. The radial dependence of b(r) for rotating stars has also these two features: the 
metric potentials of the neutron star with RMF stiff type EOS and the Q and CFL type 
quark stars give local minima for b(r), whereas the neutron star with STOS 0.5 type EOS is 
not compact enough to have a photon radius (see Table H]) . We note that the photon radii 
of the rotating stellar objects are located at higher radii than 3M, which is opposite to the 
behavior of the rotating black holes, where the photon radii approach M as the black hole 
spins up. 

The deposition rate of the total energy of the e + e~ pairs generated in the neutrino 
annihilation is usually characterized by the integral of q over the spatial proper volume of 
the neutrino stream. The integral of q in the radial direction, measuring the total amount 
of energy converted form neutrinos to electron-positron pairs at all radii R, is then given by 

r2n /*7r/2 roc 

Q{R) = 2 q(r, R, Q)^g„g e6 g <t><i) drdQd(\). (33) 

For a spherically symmetric geometry the integration over the angular coordinates <p and 
9 is a straightforward computation. This is not the case in the axially symmetric case, where 
the ^-dependence of the metric as well as the nature of the null trajectories is much more 
complicated. In order to avoid this problem, we restrict the study of the electron-positron 
energy deposition rate to neutrino pairs moving in the equatorial plane only. Instead of the 
quantity Q defined by Eq. ( |33i) . we consider its derivative with respect to in the equatorial 
plane, 

dQ f 2w f°° 

poo 

= Air q{r,R,Q)jg„ge e g H dr . (34) 
JR 9=tv/2 



dQ 
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Figure 3. The electron-positron energy deposition rate (dQ/df)) / (^dQpf/dd^j |^ in the equatorial plane around rotating 
neutron and quark stars, with the same total mass M = 1.8Mq, and same rotational velocity CI = 5 X 10 3 s — 1 . 



Q, and there: 



bre its ^-derivative, is neit her an observable, nor a Lorentz-invariant quan- 



tity (author?) fjSalmonson fc Wilson 19991 ). but Q can be used to measure the total amount 



of local energy deposited via the neutrino pair annihilation outside the neutrino- sphere. The 
derivative dQ/d9, evaluated at 9 = 71/ 2 allows one to describe the energy deposition rate 
only in the equatorial plane, but its value normalized to the Newtonian case can be applied 
to compare the energy deposition rates in different gravitational potentials. Since the quan- 
tity given by Eq. (jMj) has a simple form in the Newtonian case, the ratio of the general 
relativistic case to the Newtonian model is given by 



dQ/d9 



Ir q(r, R, 9)^g rr (r, 9)g ee (r } 9)g <M> (r, 9)dr 
e=n/2 Ir QN(r, R)r 2 dr 



=7r/2 



(35) 



dQ N /d9 

Here q(r, R, 9) is the deposition rate calculated by taking into account the general rela- 
tivistic effects, whereas ^(r, R) is the deposition rate for the simple Newtonian case, without 
taking into account the bending of the neutrino path and the redshift in the neutrino tem- 
perature. We will use the ratio given by Eq. ( 1351) to compare the electron-positron energy 
deposition rates in the space-times of neutron and quark stars, with different equations of 
state. 

In Fig. [3] we present the ratio given by Eq. ( |35i) as a function of the neutrino-sphere 
radius, measured in curvature coordinates, for each neutron and quark star model previously 
described. The total mass M of the central objects is set to 1.8M Q , and the rotational 
frequency Q of the stars is about 5 x 10 3 s _1 . For this configuration the physical parameters 
of the stars are shown in Table [TJ 

From the analysis of the compactness of these stars, that is, of the ratio of the total mass 
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EOS 


DH 


RMF stiff 


STOS 


STOS 0.5 


STOS 1 


BBBAV14 


BBBParis 


APR 


Q 


CFL150 


Pc [10 15 g/cm 3 ] 


1.29 


0.57 


0.369 


0.383 


0.40 


2.15 


1.70 


1.225 


0.931 


0.71 


M [M Q ] 


1.81 


1.80 


1.85 


1.80 


1.79 


1.80 


1.80 


1.80 


1.79 


1.80 


Mo [M @ ] 


2.05 


2.00 


2.01 


1.95 


1.93 


2.08 


2.07 


2.11 


2.09 


2.09 


He [km] 


12.01 


15.79 


21.03 


22.84 


22.72 


10.57 


10.98 


10.99 


11.79 


12.36 


f}[10 3 s _1 ] 


4.99 


5.00 


4.90 


4.71 


4.45 


5.01 


5.00 


5.00 


4.79 


5.00 


n p [10 3 s- 1 ] 


11.16 


7.97 


5.37 


4.56 


4.54 


13.96 


13.19 


13.23 


11.97 


11.28 


T/W[W- 2 ] 


3.43 


8.93 


14.91 


12.25 


9.83 


2.35 


2.66 


3.41 


4.45 


5.64 


cJ/GM? 


1.15 


2.00 


2.95 


2.49 


2.15 


0.95 


1.00 


1.12 


1.28 


1.50 


7[10 45 g cm 2 ] 


2.03 


3.52 


5.30 


4.65 


4.25 


1.66 


1.76 


1.98 


2.35 


2.63 


$ 2 [10 43 g cm 2 ] 


8.54 


43.82 


106.53 


80.44 


61.22 


4.42 


5.43 


7.87 


1.30 


1.91 


h+ [km] 


6.85 


0.00 


0.00 


0.00 


0.00 


3.19 


2.69 


0.00 


0.00 


0.00 


h- [km] 


7.48 


-3.40 


7.91 


3.87 


2.27 


8.14 


7.90 


-2.11 


0.00 


0.00 


Wc/QIlO- 1 ] 


5.85 


4.52 


4.10 


4.07 


4.08 


6.67 


6.34 


5.86 


5.27 


5.00 


Vc [km] 


9.10 


12.85 


18.00 


19.97 


19.91 


7.64 


8.06 


8.05 


8.86 


9.40 


r p /r e 


0.88 


0.72 


0.54 


0.53 


0.58 


0.92 


0.91 


0.90 


0.87 


0.84 



Table 1. Physical parameters of the compact stars with total mass M S3 1.8A/q and rotational frequency Q « 5 X 10 3 s — 1 . Here 
p c is the central density, M is the gravitational mass, Mo is the rest mass, R e is the circumferential radius at the equator, Q is 
the angular velocity, Q p is the angular velocity of a particle in circular orbit at the equator, T/W is the rotational-gravitational 
energy ratio, cJ/GMq is the angular momentum, / is the moment of inertia, $2 gives the mass quadrupole moment Mi so that 
M2 = c 4 $2/(G 2 Af 3 Mg), h+ is the height from the surface of the last stable co-rotating circular orbit in the equatorial plane, 
h— is the height from surface of the last stable counter-rotating circular orbit in the equatorial plane, oj c /0 is the ratio of the 
central value of the potential u) to f2, r e is the coordinate equatorial radius, and r p /r e is the axes ratio (polar to equatorial), 
respectively. 

to the equatorial radius, we see that the 1.8M mass stars are not compact enough to have 
a photon or a neutrino- sphere. In other words, the impact parameter b given by Eq. (1181) is a 
monotonic function of the radial coordinate, without any minimum. For comparison we also 
plotted the static case with the same total mass, which is essentially the same for any type 
of the equation of state describing the properties of the dense star matter, since the exterior 
metric of the static stars is described by the Schwarzschild metri c. The energy deposition rat e 
for this case coincides with the Q/Qn versus r/M plot given by ( Salmonson fc Wilson I999I ). 



since the multiplicative factor coming from the integration of qJg^gee9H> over # is unity 
for the Schwarzschild space-time. From Fig. [3] it can be immediately seen how the energy 
deposition rate of the electron-positron pairs is enhanced by the rotation of the central 
object. It is also clear that the measure of enhancement depends on the type of the equation 
of state used to describe the dense neutron and quark stars. In the rotating case, the increase 
in the ratio of the energy deposition rates, given by Eq. ( 1351) . is the smallest, as compared 
to the static case, for the BBBAV14 and BBBParis type EOS's. The neutron stars with 
DH and APR type EOS's produce roughly the same ratios, which are somewhat greater 
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than the BBB type neutron star pair production rates. For the quark stars, described by 
the Q and the CFL type EOS's, we obtain even higher deposition rates. Although none of 
these types of central objects possesses neutrino-spheres, they are still compact enough to 
have equatorial radii less than 5M. The equatorial radii of the neutron stars with the RMF 
and STOS type EOS's are much larger than those of the previous group: 6M for the RMF 
type EOS, and 8M or even more for the STOS type stars, depending on the temperature 
of these stellar configurations. This means that although they give the highest deposition 
rates, with the STOS T = 1 Mev type star having the maximal value, the energy released by 
the neutrino pair annihilation outside the star in the equatorial plane is smaller than in the 
case of the first group. This result might be not true in the regions close to the poles of the 
star, if the axis ratios f v j r e would be considerably smaller for the second group as compared 
to the axis ratios of the first group. In the case of the STOS type EOS f p /f e is between 0.5 
and 0.6, which is indeed much smaller than 0.85-0.9, the average values of the axis ratios for 
the first group. This means the difference between the polar radii f p of the fist group, and 
the STOS type neutron star, is smaller than the difference between their equatorial radii, 
and the lower boundary of the integral of q over the radial coordinate are closer to each 
other. In this case, the integrated deposition rate can be higher for the neutron star with 
STOS type EOS than the one with the other types of equations of state for neutron stars 
and quark stars. 

However, in this framework one should be very cautious with any statement on the 
physical processes located far from the equatorial plane, since the neutrinos reaching the 
region close to the pole of the rotating stars have impact parameters rather different from 
the neutrinos moving in the equatorial plane, and the formalism applied for the latter cannot 
be extended straightforwardly to the motion outside the equatorial plane. 

Next, we consider the electron-positron energy deposition rate for a more massive group 
of stars, with the total mass M set to 2.8M . Although this value of the mass is smaller than 
the theoretical stability mass limit of 3 Solar masses for ultra-compact objects, not all types 
of the EOS's considered here do have configurations of such high masses, even in the rapidly 
rotating case. We have obtained solutions to the field equations for this mass regime only 
for the quark stars, for the RMF, and for the STOS type neutron star EOS's, respectively. 
For the other equations of state we did not even find a compatible rotation frequency in the 
high mass regime. For the stable massive configurations the angular velocity varied up to 
6 x 10 3 s _1 for the STOS type neutron star models and the CFL type quark stars, whereas 
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EOS 


RMF stiff 


STOS 1.0 


Q 


CFL 300 


p c [10 15 g/cm 3 ] 


2.10 


0.640000 


1.60000 


0.394500 


M [M ] 


2.79986 


2.8098 


2.80017 


2.80677 


Mo [M Q ] 


3.28120 


3.19101 


3.34023 


3.28947 


i? e [km] 


14.8969 


19.6359 


15.4159 


18.0389 


f7[10 3 s- 1 ] 


9.78914 


5.90127 


10.0016 


5.96099 


n p [103s- 1 ] 


10.4349 


6.99013 


10.4087 


8.35765 


T/WIW- 1 ] 


1.30916 


1.09688 


2.00292 


1.50674 


cJ/GM^ 


5.49201 


5.05616 


6.64450 


6.02253 


7[10 45 g cm 2 ] 


4.93060 


7.52988 


5.83855 


8.87919 


* 2 [10 44 g cm 2 ] 


5.68507 


8.48428 


9.63523 


13.4904 


/i+ [km] 


1.93604 


0.306693 




0.00000 


h- [km] 


19.6015 


15.1710 


21.4442 




wc/nficr 1 ] 


8.23642 


6.2819 


7.57682 


5.75394 


r e [km] 


9.99343 


14.9859 


10.2897 


13.1041 


r p /r e 


0.590234 


0.640000 


0.502344 


0.620000 



Table 2. Physical parameters of the compact stars with total mass M « 2.8Af0 and rotational frequency f2 « 6 X 10 3 s 1 and 
10 4 s" 1 . 



EOS 


RMF stiff 


STOS 1.0 


Q 


CFL 300 


M/R e 


0.2775 


0.2113 


0.2682 


0.2298 


R s /M 


3.63 


4.73 


3.73 


4.35 


Re [km] 


14.89 


19.64 


15.42 


18.04 


R P h/M 


2.14 




2.88 


4.24 


Rph [km] 


8.82 




11.87 


17.56 



Table 3. The compactness, the equatorial and the (equatorial) photon radii of the compact stars with total mass M «s 2.8Mq 
and rotational frequency fi»6x 10 3 s —1 and 10 4 s . 

Q reached values of around 9 x 10 3 — 10 4 s _1 for the neutron stars with the RMF type EOS, 
and for Q type quark stars. The models we considered for the study of the electron-positron 
energy deposition rate have a rotational frequency of 6 x 10 3 s -1 , for the first two models, 
and a frequency of 10 4 s _1 for the second ones, as shown in Table [2J These configurations, 
except the model with the RMF type EOS, are already compact enough to have photon or 
neutrino radii in the equatorial plane. The values of the neutrino radii together with the 
equatorial radii of the stars are given in Table HI The radial profile of the ratio in Eq. (|35|) 
for these configurations, and of the static case, is presented in Fig. HI 

The ratios {^dQ/dO^ / {dQ^/dO^j ^ ^ obtained for the CFL quark stars are somewhat 
higher than those obtained for the STOS type EOS's, and this ratio is essentially 45-50% 
higher at the radius R = 5M, as compared to the same quantity for the static stars. As 
for the central objects with RMF and Q type EOS's, the rotating neutron star produces a 
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Figure 4. The electron-positron energy deposition rate {dQ/dd^ / (dQw/dO^ in the equatorial plane around rotating 

neutron and quark stars, with the same total mass M = 2.8M , and rotational velocities Q = 6 X 10 3 s" 1 (for STOS and CFL 
type EOS's), and f2 = 10 4 s" 1 (for RMF and Q type EOS's). 

deposition rate about 35% higher than in the static case. For the quark star at the radius 
of R = 5M, the increase in the deposition rate can reach 50% more than in the static case. 
In the region R w 5M, close to the surface of the rotating neutron and quarks stars, and 
which is relevant for both the neutron star collapse and the supernova explosion, there is a 
considerable increase in the energy deposition rate of the electron-positron pair creation, as 
compared to the Newtonian case. 

We have compared not only the neutrino-antineutrino annihilation for different types of 
central objects, but we have also studied the effects of the different rotational velocities of 
the same central object on the energy deposition rate. We have chosen the APR EOS for 
neutron stars, and the Q EOS for quark stars, to produce sequences of rotating stars with 
different angular velocities, but with a fixed total mass of 1.8M . The rotational frequencies 
Q are shown in Table HI together with the other physical parameters of these configurations. 



In Fig. [5]we present the ratio (dQ/d9) / idQ^/dO 



= 7r/2 



calculated for configurations 



with increasing angular velocities. With increasing VL the deposition rate is also increasing for 
both types of central objects. The ratio (dQ/dQj / {dQ^/d9j ^ ^ reaches 15% more for the 
quark star than its value in the static the angular frequency approaches the Keplerian 

limit. Since the neutron star rotating with the angular velocity of 8 x 10 3 s _1 becomes less 
compact, its equatorial radius is somewhat less than 7M, and the deposition rate cannot be 
integrated to lower radii. This prevents the star from accumulating much energy released 
in the form of electron-positron pairs, even if the ratio given by Eq. (1351) is very high in 
the equatorial plane for R > 7M, as compared to the static case. However, these results 
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Table 4. Physical parameters of the neutron stars with APR type EOS and Q type quark stars, with a total mass of M «3 1.8Mg, 
and with different rotational frequencies. 




R/M R/M 

Figure 5. The electron-positron energy deposition rate [dQ/dd} / {^dQ^/dO^ | ^ in the equatorial plane of a rotating 

neutron star, with the APR EOS (left hand side), and for a quark star with EOS Q (right hand side), the stars having the 
same total mass M = 1.8Mq, but different rotational velocities. 

show that the rotating ultra-compact objects can convert some of their rotational energy 
into the energy of the e + e~ pairs. The faster the central objects rotates, the more energy 
will be released in the neutrino-antineutrino pair annihilation process. Then we expect that 
ultra-compact stellar object, especially quark stars, rotating close to the Keplerian limit, can 
deposit a large amount of energy in their equatorial plane due to the neutrino-antineutrino 
annihilation process. 

Similarly to the case of the energy deposition, one can define the integrated momentum 
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Figure 6. The electron-positron momentum deposition rate (dP/dd^ / (dP^/dO^ |^ in the equatorial plane around ro- 
tating neutron and quark stars, with the same total mass M = 1.8Mg, and same rotational velocity f2 = 5 X 10 3 s — 1 . 



deposition rate as 

r2ir /■tt/2 roo 



rZTT /"7T / A pOO 

P(R) = 2 J y J r p{r,R,6)^g^g^drd9d0, (36) 
where p is given by Eq. (pD). In the equatorial plane we introduce the ratio 
dP/dO 



lRp{r, R, 6)^g rr {r, 0)gee{r, 9)g H> {r, 0)dr 



= 7r/2 



0=7r/2 J^p N (r,R)r 2 dr 



(37) 



dP N /d6 

of the momentum deposition rate in the general relativistic case and of momentum deposition 
ratio P/v corresponding to the Newtonian case. In Fig. (El) we present the plots of this 
quantity as a function of R/M for neutron and quark stars with 1.8M , and with a rotation 
velocity of 5 x 10 3 s~ 1 , respectively. This figure is a counterpart of Fig. (j3J), displaying the 
energy deposition rate for the same configuration of the stellar objects. By comparing the two 
figures we see that they exhibit the same characteristics for each stars, i. e., the momentum 
deposition rate is roughly proportional to the the energy deposition rate. Therefore the RMF 
and STOS EOS type stars produce the highest ratios for the momentum deposition rate. 
Quark stars have a smaller momentum deposition rate, and we obtain the smallest values 
for the neutron stars with DH, APR and BBB type EOSs, respectively. 



6 DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered the general relativistic effects on the energy deposi- 
tion rate from the neutrino-antineutrino annihilation process near rapidly rotating neutron 
and quark stars. The energy deposition rate has been obtained numerically for several equa- 
tions of state of the neutron matter, and for two types of quark stars, respectively. All the 
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general relativistic correction factors, related to this process, can be obtained from the met- 
ric of the central compact object. Due to the differences in the space-time structure, the 
quark stars present some important differences with respect to the energy deposition rate, 
as compared to the neutron stars. As a general result we have found that rotation always 
enhances the energy deposition rate. 

A possible astrophysical application of our results could be in the explanation of the 
physical processes that lead to the formation of the Gamma- Ray Bursts (GRBs). The so- 
called fireball model can basically explain the observational facts well, and thus it is strongly 
favored, and widely accepted today (for recent reviews on GRBs see Meszaros (2006) and 
Zhang (2007), respectively). In this model, the central engine gives birth to some energetic 
ejecta intermittently, like a geyser, producing a series of ultra-relativistic shells. The shells 
collide with each other and produce strong internal shocks. The highly variable 7-ray emis- 
sion in the main burst phase of GRBs should be produced by these internal shocks. After 
the main burst phase, the shells merge into one main shell and continue to expand outward. 
These merged shells, which moves relativistically, interacts with the surrounding medium to 
form an external shock. This shock will give a good fraction of its energy to the swept-up 
electrons and accelerate them to relativistic velocity. Similarly, a fraction of shock energy 
will go to the magnetic field. These shocked relativistic electrons move in the magnetic field 
and emit synchrotron radiation to produce a broadband electromagnetic emission called 
"afterglows". According to the Swift observations(see Zhang (2007) for a review), the decay 
of the X-ray afterglow can be classified roughly into four stages. Stage I (up to ~ 100 s) 
indicates the early fast decay phase, usually with flux density decay as t~ 3 to t~ 5 ; Stage II 
(up to ~ 1000 s) indicates the subsequent shallow decay phase, with t~ 2 to t -0 ' 8 ; Stage 
III (up to ~ 10 4 s) is the late normal decay phase, t _1 to t _L3 ; Stage IV (beyond 10 4 s) is 
the late fast decay phase, t~ 2 to t -2 5 . The first stage is explained as the tail decay of the 
prompt emission, and the third and the fourth stages can be explained very well in terms of 
the standard fireball model. Although there are many different explanations of the second 
stage, the "shallow decay phase", one of the most popular explanations is the late time 
continuous energy injection. Obviously, the energy deposition in the equatorial plane cannot 
be responsible to the prompt gamma-ray emission, which requir es a fast and a very narrow 
beam jet. In fact, the high energy photons detected by Fermi (IMeegan et al. 20091 ) imply 
that the Lorentz factor of fast jet is over 500. On the other hand, most deposition energy 
should come from the equatorial plane. We suggest that the energy deposition along the axis 
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will form a narrow cone jet, which encounters much less material, and hence it can maintain 
a very large Lorentz factor. On the other hand, the energy deposition in the equatorial plane 
is higher, but it also encounters much more material, and hence it moves slowly When the 
fast jet slows down due to the interstellar medium, this slow but more energetic wind from 
the equatorial plane can inject additional energy into the fast jet, and result in the shallow 
decay phase. The detailed calculations will be carried out in our future study. 

The neutrino annihilation and the electron-positron pair production also plays an essen- 
tial role in the astrophysical processes associated with the phase transitions that could 
take place inside neutron stars. For example, the sudden phase transition from normal 
nuclear matter to a mixed phase of quark and nuclear matter induces temperature and 
density oscillations at the neutrinosphere. Consequently, e xtremely intense, p ulsating neu- 
trino/antineutrino and leptonic pair fluxes will be emitted (jCheng et al. 20091 ). During this 
stage several mass ejecta can be ejected from the stellar surface by the neutrinos and antineu- 
trinos. These ejecta can be further accelerated to relativistic speeds by the electron/positron 
pairs, created by the neutrino and antineutrino annihilation outside the stellar surface. In 
order to produce the Gamma- Ray Bursts, a high neutrino emission rate is necessary. On the 
other hand, it is important to note that electron-positron pairs can deposit energy much more 
efficiently than the neutrinos, and the dominant ene rgy deposition process is the neutrino- 
antineutrino annihilation process (jCheng et al. 20091 ) . In fact most pairs are created outside 
the star. A large fraction of the neutrino energy, will be absorbed by the matter very near 
the stellar surface. When this amount of energy exceeds the gravitational binding energy, 
some mass near the stellar surface will be ejected, and this mass will be further acceler- 
ated by absorbing pairs created from the neutrino and antineutrino annihilation processes 
outside the star. T his process may be a possible mechanism for short Gamma-Ray Bursts 
( Cheng et al. 2009^ 

The neutrino emission rate is also strongly dependent on the temperature. In the stan- 
dard models of GRBs it is assumed that the central object is surrounded by a degenerate 
accretion di sk, which allo ws super-Eddington accretion rates, of the order of one solar mass 
per second (jZhang 20071 ). If the central compact object is a neutron or a quark star, such 
super-Eddington accretion rates can maintain the compact object at very high tempera- 
tures, and thus allowing very high neutrino luminosities, as well as a high rate of electron- 
positron pair production. The neutrino temperature can be estimated by assuming that 
the accretion power r]Mc 2 , where r\ is the efficiency of the energy conversion and M is 
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the accretion rate, is equal to the radiation power AiiR 2 c>T^ which gives the temperature 
as T v = (r]Mc 2 /AixR 2 a^ ^ . By taking r\ = 0.1, R = 10 6 cm, and an accretion rate of 
M = 1Mq/10 s, we obtain T v = 7.14 x 10 10 K, a temperature which is of the order of MeV. 
Therefore, if the accretion disk is fed at a high rate, like, for example, by the fallback material 
after a supernova explosion, a high neutrino-antineutrino emission rate can be maintained, 
and this could explain some of the basic properties of the GRB phenomenon. 
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